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We investigate the finite size effect on masses of excited baryons in quenched lattice QCD simula- 
tion. For this purpose, we perform numerical simulations at three different lattice sizes. La ~ 1 .6, 
2.2 and 3.2 fm. The gauge configurations are generated with the single plaquette gauge action 
at j3 = 6/§^ = 6.2, and the quark propagator are computed with the Wilson fermion action. To 
access to two parity states of the nucleon and four different spin-parity states of the A baryon, the 
appropriate spin/parity projection are carried out. We find that the spatial lattice size is required 
to be as large as 3 fm to remove the finite size effect on excited baryons even in the heavy quark 
region (Mj^/Mp ~ 0.82, 0.87), where that on the nucleon is negligible. On our largest lattice 
{La ~ 3.2 fm), all mass spectra of the 7^ = 1 /2^ nucleons and the 7^ = 1 /2^,3/2* A baryons 
are roughly consistent with experimental values after naive chiral extrapolation. 
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Recently, many lattice studies on excited baryon spectroscopy have been made within the 
quenched approximation, and gave results in gross agreement with experiments [[I), ^, ^, ||, |^ 0]. 
To achieve the higher precision study, however, it is necessary to control several systematic errors, 
which may arise from the quenched approximation, chiral extrapolation, nonzero lattice spacing 
and finite size effect. In this study, we especially focus on the finite size effect, which has been not 
much paid attention so far. 

From the phenomenological point of view, the finite size effect may be interpreted by the 
"wave function" squeezed in the finite volume [||]. Accordingly, to remove this effect, the wave 
function should be completely fitted in a finite lattice box. As reported in Ref. [^, the spatial size 
of lattice is required to be La > 2.5 fm even for the ground-state nucleon in quenched lattice QCD. 
It is naively expected that the much larger size is required for the excited state than the ground 
state. However, most of the past studies are performed with spatial lattice sizes La ~ 1.6 — 2.2 fm 
due to large computational costs ' . 

To estimate the finite size effect, we perform simulations with three different lattice sizes. 
La ~ 1.6, 2.2 and 3.2 fm. We examine this effect on masses of the ground and excited baryons, 
the 7^ = 1 /2^ nucleons and the = 1 /2^ ,3/2^ A baryons. In general, the local baryon operator 
couples to both the positive- and negative-parity state. The appropriate parity projection is desired 
since the energy level of the negative-parity baryon is higher than that of the ground state. The 
spin projection is required for a precise treatment of the A baryon spectra. We briefly explain the 
projection technique that we utilized in this study as below. 

The local baryon operator &b{x) = £abc{qi{^)CTqi,{x))T' qc{x) couples to both the positive- 
and negative-parity state regardless of the intrinsic parity of the operator [|l|]. Thus, the two-point 
function Gb(0 = "Li{0\T {0'B{x,t)^B(O,O)]\O) has the following asymptotic form [[IHl^, 

Gs{t) =A,i±^f(^.-^^l'l _AJ-^^^e-"-\^\ (1) 

at large Euclidean time t. Here A-t and M± denote the amplitude and mass of the lowest-lying state 
in each parity channel, respectively. 

On the lattice with finite extent T in the time direction, the two-point function receives the 
reflection from the time direction boundary. Under the periodic/anti-periodic boundary condition, 
Eq. dl]) should be expressed as 



^p.b.c/a.p.b.c.^^^ 



+ 2 2 



(2) 



where only the primal reflections are included since the higher orders of wrap-round effect are 
negligible in comparison with the first wrap-round effect. Clearly, unwanted-parity contamina- 
tions originated from the first wrap-round effect cannot be removed by the simple projection 
Tr I^.qP^ '^ I^ 'p-^ '^- (^i'^ Xo extract the state with desired parity, we employ the linear combi- 
nation of two correlation functions with periodic and anti-periodic boundary conditions, GB{t) = 



'in Ref.[^, a large lattice size {La ~ 3.2 fm) for the simulation is achieved with a relatively coarse lattice spacing. 
We however recall that masses of excited baryons are beyond their cutoff scale (a^' ~ 1 GeV). 
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l^p.b.c. + Gg'' "' '^' (f)}/2. As the results, the free-boundary situation is realized, so the simple 
projection becomes efficient. In the numerical simulation, we adopt a linear combination of the 
quark propagators. As shown in Ref. |11], the linear combination in the quark level automatically 
realizes the linear combination in the hadronic level. We perform this parity projection to apply 
to baryon two-point correlators. We access to the two different parity states of the nucleon, which 
correspond to the positive-parity state (A'^) and the negative-parity state (A^*). 

We also treat the A baryon in this study. The simplest composite operator for the A baryon 
(yf{x) = EabA^a^^^yi^b{^\)Uc{x) has the structure of the Rarita-Schwinger spinor so that the two- 
point function constructed from this operator coupled to both the spin-3/2 and spin- 1/2 state [p^]. 
Here / is the spatial Lorentz index. The two-point function Gfy(0 = Lv(0|7'{^f (^>0^f (0>0)}|0) 
is expressed by the orthogonal sum of spin-3/2 and spin- 1/2 components: 

Gf;(0 = (5,7 - \yiy^ ^3/2(0 + \ymG\i2^t) o) 

with appropriate spin projection operators [|I3|]. Therefore, we can access the excited A state such 
as the spin- 1/2 A state by the spin projected correlators 

G^/2(0 = |Gt(0-^lT^T.Gt-(0, (4) 
z z ^ 

G\|^{t) = Y^yakG^kSf)- (5) 

k 

We perform the appropriate parity projection to apply to G'^i^it) and G'^i^(t) as well. We finally 
access to the four different spin-parity states of the A baryon, which correspond to the spin-3/2 
positive-parity state (A3/2), the spin-3/2 negative-parity state (A3^2)' the spin-1/2 positive-parity 
state (A1/2) and the spin-1/2 negative-parity state {t^yj-^- The A3/2 is the ground state of the A. 

The gauge configurations have been generated in quenched approximation with the single 
plaquette gauge action at j8 = = 6.2, where the cut-off scale (~ 3 GeV) was definitely higher 
than mass scale of all calculated baryons (~ 1 — 2 GeV). To examine the finite size effect, we 
performed simulations on three different lattice sizes, x T = 24^ x 48 (A'conf = 320), 32^ x 
48 (A'^conf = 240) and 48^ x 48 (A'^conf = 210). To compute the quark propagator, we utilized the 
standard Wilson fermion action, of which computational costs are relatively cheep. The hopping 
parameter K- = 0.1520, 0.1506, 0.1497, 0.1489, 0.1480 cover the range M;r/Mp ~ 0.66 - 0.96. 
To perform the appropriate parity projection, we adopt a procedure to take an average of two quark 
propagators, which are subjected to periodic and anti-periodic boundary condition in time. We 
calculated the point-to-point quark propagator and constructed the baryon two-point functions by 
the conventional interpolating operators, eabc{ulCY5dh)uc for the nucleon, and £abciulCYiUi,)uc for 
the A baryon. 

In Fig.|l|, the lattice-size dependence on the mass of each baryon (N, N*, A3/2, A3/2) shown 
for two hopping parameters (k = 0.1506 and 0.1497), which correspond to the relatively heavier 
quark masses. The quoted errors represent only statistical errors. Horizontal dashed lines represent 
the values in the infinite volume limit together with their one standard deviation. Those values are 
evaluated with the phenomenological power-low formula, aMi = aMco + cL^^. 

For the ground-state of the nucleon, we do not observe any serious finite size effect even on 
our smallest lattice. On the other hands, we find that the mass of the negative-parity nucleon suffers 
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Figure 1: Masses of the = 1/2^ nucleons (left figures) and the = 3/2* A baryons (right figures) as 
function of spatial lattice size in lattice units for two hopping parameters (k = 0. 1497 and 0. 1506). 
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Figure 2: A' — A3 /2 hyperfine splittings (left figures) and A^* — Ay 2 hyperfine splittings (right figures) as 
function of spatial lattice size in lattice units for two hopping parameters (jc = 0. 1497 and 0. 1506). 



from the large finite size effect. The observed tendency of the size effect is that the mass decreases 
as the lattice size increases. This result seems to be consistent with the picture where the squeezed 
"wave function" may increase the total energy of the three-quark system and the "wave function" 
of the excited state has the larger spatial extent than that of the ground state. 

For the A3/2 baryon, however, we find that the finite size effect has a different pattern: Ma, , 
becomes large as the lattice size increases. On the other hands, the Ay^ has the same pattern of 
the finite size effect observed in the A^* spectrum, while the finite size correction of the is 
relatively milder than that of the A^*. We confirm the similar feature in the case of the A1/2 baryon 
as well. In all cases, to keep the size effect at a few percent levels, the spatial size of lattice is 
required to be at least 3 fm. 

The peculiar behavior on the A3/2 may originate from a hyperfine interaction. Our simulations 
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are carried out in the relatively heavy quark region, where the size effect of the nucleon is almost 
negligible. Thus, the finite size effect on the A3 /2 can attribute to the effect of the A'^ — A3 /2 hyper- 
fine mass splitting. In Fig. ^ the lattice-size dependence of the hyperfine mass splitting in each 
parity channel is shown. In the negative-parity N* — A*y2 hyperfine splitting, we observe the same 
pattern as the N — At,^2 hyperfine splitting. The same feature is also observed in the charmonium 



spectrum [14]. The hyperfine splitting between rjc and J/xj/ diminishes as the spatial lattice size 
decreases if La < 1.3 fm. It seems that it is difficult for one-gluon exchange as a spin-spin compo- 
nent of the Fermi-Breit type interaction to account for this particular pattern of the finite size effect 



on the hyperfine interaction [11] 



Next, we extrapolate the baryon masses calculated on our largest lattice {La ~ 3.2 fm) toward 
the chiral limit. Our quark masses are too heavy to employ guidance by the chiral perturbation 
theory, so the empirical curve fit formula is alternatively utilized: 

(uMb)^ = dQ + diiaMnf. (6) 

As reported in some Refs. [^, [Tl|], this formula gives a better fit to the data in the heavy quark 
region than the linear fit. 

In Table |TJ the baryon masses in the chiral limit are listed in the physical unit. Two types of 
input (ro input and Mp input) are taken to reveal the dependence of the choice of input to set a scale. 
Since the systematic errors stemming from this choice exceed the statistical errors of the nucleon 
or the A3/2, we alternatively quote various mass ratios also. 



N 



N 



1.28(4) 


(Expt. : 


-1.31), 


1.61(10) 


(Expt. : 


~ 1.63), 


1.28(7) 


(Expt. : 


~1.11), 


1.61(7) 


(Expt. : 


~ 1.38). 



(7) 



The mass ratios show a good agreement with the experimental values within statistical errors except 
for the case where Ma*^^ is included. The values extrapolated from large quark mass should not 
be taken seriously. However, the level ordering in A spectra, Mt^^^^ < Ma*^^ ^ ^^'3/2 ^ ^^1/2' 
well-reproduced in comparison to corresponding A states, which are all ranked as four stars on the 
Particle Data Table [p3|]. In addition, it is worth mentioning that a signal A( 1750) (7 = 3/2 and = 
1/2+), which is the weakly established state (one star) |[T^], cannot be seen in our data. 

In this study, we have investigated the finite size effect on masses of the ground and excited 
baryons, the = 1/2^ nucleons and the 7^ = 1/2^ , 3/2^ A baryons. For this purpose, we perform 
numerical simulations at three different lattice sizes. La ~ 1.6, 2.2 and 3.2 fm. We have found the 
sizable finite size effect on masses of all A states and the negative-parity nucleon even in the heavy 
quark region (Mj^/Mp ~ 0.82, 0.87), where that on the nucleon is negligible. To keep this effect 
at a few percent levels, the spatial size of lattice is required to be at least 3 fm. We found a 
peculiar lattice-size dependence on the mass of the positive-parity A3/2 state: the mass of the A3/2 
state becomes small as the lattice size decreases. This behavior seems to originate from the finite 
size effect on the hyperfine interaction. We also extrapolate the baryon masses toward the chiral 
limit with the empirical curve-fit formula. All mass spectra of the = 1/2* nucleons and the 
7^ = 1/2^,3/2* A baryons are roughly consistent with experimental values. 
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Table 1: The third and forth column list results of all measured baryon masses in GeV units, which are set 
by two different inputs, input and Mp input. The fifth and sixth list experimental values of the correspond- 
ing baryon and its status in the Particle Data Table JT^. The possible assignments of the SU{6) (8)0(3) 
supermultiplet are also embedded into the final column. 



baryon 




our results [GeV] 
(ro input) 


{Mp input) 


physical state (status) 


5^7(6)0 0(3) 
classification 




V2' 2 ! 


0.968(23) 


1.031(25) 


A^(939) **** 


[56,0+] 




V2' 2 > 


1.555(97) 


1.658(102) 


A^(1535) Sn **** 


[70,1-] 


A3/2 


(3 3 + ^ 
V2' 2 1 


1.236(29) 


1.320(31) 


A(1232) **** 


[56,0+] 


A* 
'^3/2 


(I l-\ 
V2' 2 ) 


1.985(72) 


2.114(75) 


A(1700) D33 **** 


[70,1-] 


Al/2 


(3 i+^ 
V2' 2 ) 


2.325(112) 


2.478(118) 


A(1910) ^31 **** 


[56,2+] or [70,0+] 


A* 
'^1/2 


ll l-\ 
V2' 2 ! 


1.866(153) 


1.987(165) 


A(1620) 531 **** 


[70,1-] 
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